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Abstract 

We report the results of a study of B ± — > K ±r q c and B ± — > K ±r q c (2S) decays fol- 
lowed by r] c and n c (2S) decays to (KsKit) . The results are obtained from a data 
sample containing 535 million i?i?-meson pairs collected by the Belle experiment 
at the KEKB e + e~ collider. We measure the products of the branching fractions 
B(B ± -> K ± r ?c )B(r ?c -> KsK^) = (26.7±lA(stat)tli(syst)±4.9{model)) x 1(T 6 
and B{B ± -> K ± 7 lc (2S))B(7] c (2S) -»• KgK^) = {SA+f^stat+model^oiisyst^x 
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10 . Interference with the non-resonant component leads to significant model un- 
certainty in the measurement of these product branching fractions. Our analysis 
accounts for this interference and allows the model uncertainty to be reduced. 
We also obtain the following charmonia masses and widths: M(j) c ) = (2985.4 ± 
1.5(stat)+^(syst))MeV/c 2 ,r(r/ c ) = (35.1±3.1(stat)t{i(syst)) MeV/c 2 , M(r? c (25)) = 
(3636.11^ (stat+model)±% 7 (syst)) MeV/c 2 , T(r] c (2S)) = (6.6t 8 5 t(stat+model)tH(syst)) 
MeV/c 2 . ' 

Key words: B decay, etac 

PACS: 13.25.Gv, 13.25.Hw, 14.40.Pq 



1 Introduction 

Charmonium states consist of a heavy charm- ant icharm quark pair, which 
allows the prediction of some of the parameters of these states using non- 
relativistic and relativistic potential models [1], lattice QCD [2], non-relativistic 
effective field theory (NRQCD) [3], and sum rules [4] (see the recent review 
in [5]). The comparison of these predictions with experimental results provides 
an opportunity to tune the parameters of theoretical models and, therefore, 
improve the accuracy of other values predicted by these models. We have to 
measure the charmonium masses, widths, and product branching fractions 
with enough accuracy to compare them with theoretical predictions. Param- 
eters of (cc) states such as the r] c and i] c (2S) mesons have been studied in 
various experiments using a variety of decay channels [6]. Tables 1 and 2, 
which list a selection of most precise mass and width determinations, show 
that there is quite a large spread of measured masses and widths of the r\ c and, 
especially, i] c (2S) mesons resulting in large scale factors for the world average 
values [6] . Moreover, our knowledge of hadronic decays of these charmonia is 
rather poor. 

In the Belle and BaBar B factory experiments, charmonia are produced in 
various ways: from fragmentation in electron-positron annihilation, from two- 
photon processes, and in B decays. The advantages of B ± — >■ K^cc decay 
are the relatively large reconstruction efficiency, small background, and the 
fixed quantum numbers (J p = 0~) of the initial state. Here we consider the 
following decays of charged B mesons: 

B ± ->■ -> ^(KsKitf, 

B ± ->■ X ± ?7 C (25) ->■ ^(KsKnf. 

A 492 fb" 1 data sample provides an opportunity to determine the correspond- 
ing products of the branching fractions as well as masses and widths of the rj c 
and 7) C (2S) mesons. 
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Table 1 

Previously measured r/ c parameters. 



Experiment 


Process 


Mass, MeV/c 2 


Width, MeV/c 2 


BaBar [7] 
Belle [8] 
BaBar [9] 
Belle [10] 
CLEO [11] 
BES [12] 
E835 [13] 


77 ->■ Vc -> K^K^ir^ 
77 ~~ ^ Vc — > KsK^ic^ 
B -> ^ifW -»• KKttK^ 
77 - ^ — 7-liadrons 
77 -> % ->■ -K'°-ftT ± vr =F 
J/^ ->• 7r/ c 
77 


2982.2 ±0.4 ±1.6 

2981.4 ±0.5 ±0.4 
2985.8 ±1.5 ±3.1 
2986.1 ± 1.0 ± 2.5 
2981.8 ± 1.3 ± 1.5 

2977.5 ± 1.0 ± 1.2 
2984.1 ±2.1 ±1.0 


31.7 ±1.2 ±0.8 
36.6 ±1.5 ±2.0 

36.3l|£ ± 4.4 
28.1 ±3.2 ±2.2 

24.8 ± 3.4 ± 3.5 
17.0 ±3.7 ±7.4 

20.4+^ ± 2.0 



Table 2 

Previously measured rj c (2S) parameters. 



Experiment 


Process 


Mass, MeV/c 2 


Width, MeV/c 2 


Belle [8] 
Belle [14] 
BaBar [15] 
CLEO [11] 
BaBar [16] 
Belle [17] 


77 -> r] c (2S) -»• KgK^it^ 
e + e~~ — >■ J/i/jcc 
e + e~ — >■ J/tpcc 
77 -> r/ c (2S") -»• ifgif ±7^ 
77 -»• »/ c (25) -»• KKtt 
B -»• KK s K±iF 


3633.7 ±2.3 ±1.9 
3626 ± 5 ± 6 

3645.0 ± 5.5+^| 
3642.9 ±3.1 ± 1.5 

3630.8 ± 3.4 ± 1.0 
3654 ± 6 ± 8 


19.1 ±6.9 ±6.0 

22 ± 14 
6.3 ± 12.4 ±4.0 
17.0 ±8.3 ±2.5 
< 55 



At all stages of this analysis we consistently take into account the interference 
between the B ± — > K ± r/ C and B ± — > K^rf^S) decays and the decay B ± — > 
K ± (KsKn)°, which has the same final state but no intermediate charmonium 
particle. 



2 Event selection 



The results are based on a data sample that contains 535 x 10 6 BB pairs, 
collected with the Belle detector at the KEKB asymmetric-energy e + e~ col- 
lider [18] operating at the T(AS) resonance. 

The Belle detector [19] is a large-solid-angle magnetic spectrometer that con- 
sists of a silicon vertex detector (SVD), a 50-layer central drift chamber (CDC) 
for charged particle tracking and specific ionization measurement (dE/dx), an 
array of aerogel threshold Cherenkov counters (ACC), time-of-flight scintilla- 
tion counters (TOF), and an array of 8736 CsI(Tl) crystals for electromag- 
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netic calorimetry (ECL) located inside a superconducting solenoid coil that 
provides a 1.5 T magnetic field. An iron flux return located outside the coil 
is instrumented to detect K\ mesons and identify muons (KLM). We use a 
GEANT-based Monte Carlo (MC) simulation to model the response of the 
detector and determine its acceptance [20]. 

Pions and kaons are separated by combining the responses of the ACC and 
the TOF with dE/dx measurements in the CDC to form a likelihood C(h) 
where h = n or K. Charged particles are identified as pions or kaons using 
the likelihood ratio 1Z: 



Charged tracks are selected with requirements based on the x 2 °f the track 
fits and the impact parameters relative to the interaction point. We require 
that the polar angle of each track be in the angular range 18° - 152° and that 
the track momentum perpendicular to the positron beamline be greater than 
100 MeV/c. 

Charged kaon candidates are identified by the requirement 7Z{K) > 0.6, which 
has an efficiency of 90% and a pion misidentification probability of 3 - 10% 
depending on momentum. For pion candidates we require TZ(tt) > 0.2. K s 
candidates are reconstructed via the n + 7i~ mode. We apply the following cut 
on the 7T+7T- invariant mass: 0.489 GeV/c 2 < M(it + n-) < 0.505 GeV/c 2 . The 
flight length of the Ks is required to lie within the interval [0.1 : 20] mm. The 
condition on the Ks direction angle f is cos(ip) > 0.95. 

B meson candidates are identified by their center-of-mass (cm.) energy differ- 
ence AE = (Y,i Ei)—E h , and the beam-constrained mass M bc = \J E^ — (Y^iPi) 2 , 
where E\> = yfs/2 is the beam energy in the T(4S) cm. frame, and and 
Ei are the cm. three- momenta and energies, respectively, of the B meson 
candidate decay products. The signal region is defined as: \AE\ < 0.03 GeV, 
5.273 GeV/c 2 < M bc < 5.285 GeV/c 2 . The AE sideband region is defined as 
\\AE\ -0.06(GeV)| < 0.03 GeV. 

To suppress the large continuum background (e + e~ — > qq, where q — u,d, s, c), 
topological variables are used. Since the produced B mesons are nearly at 
rest in the cm. frame, the signal tends to be isotropic while continuum qq 
background tends to have a two-jet structure. We use the angle between the 
thrust axis of the B candidate and that of the rest of the event (©thrust) 
to discriminate between these two cases. The distribution of I COS ©thrust IS 
strongly peaked near | cos ©thrust | = 1 for qq events and is nearly uniform for 
T(4S') — > BB events. We require | cos ©thrust | < 0.8. 
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If there is more than one combination that satisfies the selection criteria, we 
select the B candidate with the minimum difference \M{Ks) — M(7r + 7r~)| 
and the minimum difference between the vertex z-coordinates of the kaon and 
pion from the charmonium decay, and the kaon from the B decay. If the final 
state includes two kaons of the same charge, as in K (K Kstt^), we must 
choose the one from the charmonium decay. For this purpose we select the 
candidate with the minimum difference \M(r) c /r) c (2S)) — M(KsK ± tt t )\. The 
mean number of multiple candidates per event is 1.6. 



3 Interference study 

The data sample can contain signal, which has the same final state 
onant decay but does not include a charmonium resonance. The contribution 
of these events is referred to as the non-resonant amplitude. Since the final 
state is the same, this amplitude interferes with the signal. However, if the 
final particles form narrow resonances such as D, D$, and (p mesons, the in- 
terference effect cancels after integration over mass. Thus, we can reject some 
background channels by applying appropriate cuts on the mass combinations 
of the final state particles. If the intermediate resonances have a substantial 
width or the B meson decays directly into the final state particles, the effect 
of interference must be taken into account. 

The non-resonant contribution can be seen as a peak in the AE distribution 
in the charmonium sideband regions shown in Figs. 1 and 2 (plots on the 
right jjj. We fit the AE distributions with the sum of a Gaussian distribution 
and a second-order polynomial function. From these fits we obtain the number 
of events in the signal region N f, s and in the sideband region N s t>, which can 
be rescaled to obtain the number of non-resonant events N non _ res . In the 7] c 
case N b s = 889 ± 37(stat) and N non - res = 87 ± ll(stat). In the r] c (2S) case 
N obs = 279 ± 29(stat) and N non _ res = 156 ± 13(stat). 

Different values of the interference phase can give significant variations in the 
number of signal events while the total number of observed events remains 
the same. Using the hypotheses of maximal constructive and destructive in- 
terference, we would obtain 410 and 1550 signal events, respectively. There- 
fore, the model uncertainty in the number of rj c signal events is rather large: 
N S ig n ai = 980 ± 570(model). It would be even larger for the T] C (2S) decay. A 
dedicated study of the interference effect allows this uncertainty to be reduced. 



1 The rj c meson mass signal region is (2.92 - 3.04) GeV/c 2 , the sideband is (2.54 
- 2.86) GeV/c 2 and (3.14 - 3.46) GeV/c 2 . The rj c (2S) meson mass signal region is 
(3.58 - 3.7) GeV/c 2 , the sideband is (3.14 - 3.46) GeV/c 2 and (3.72 - 4) GeV/c 2 . 



7 




AE, GeV AE, GeV 



Fig. 1. AE distributions in the signal (left) and sideband (right) regions for 
B ± -> K±i\ c -> ^{KsKtt) candidates. 




Fig. 2. AE distributions in the signal (left) and sideband (right) regions for 
B ± ->• K ± ? ?c (25') -»• K+iKsKir) candidates. 

In the B — > KmKgK^Tt decay (see Fig. 3) there are four particles in the final 
state, which gives 4x3 measured parameters. Taking into account the four 
constraints of energy-momentum conservation and integrating over the three 
angles that characterize the B decay (it is a pseudoscalar and there should be 
no dependence on these angles), we have 5 independent variables to describe 
the amplitude of the process. We chose the following parameters: KgKtyir 
invariant mass, two Dalitz variables for the r] c (or i] c (2S)) decay - q\ and q\ 
(for example, M(K(2)7i) 2 and M(Kstt) 2 ), the angle between the K$ and 
in the rest frame of the K^Kg'K system (8), and the angle between the planes 
of the K(i) — n and K(i) — Ks in the same system (0). 

The M(K s Ktt) distribution has four peaks corresponding to r) c , J/ip, Xd 
and r) c (2S) production (see Fig. 4). In addition to these peaks, there is a 
non-resonant signal, which interferes with the i] c (or r/ c (25)) signao. Unfortu- 
nately, the shape of the one dimensional (1-D) mass distribution alone does not 
allow the interference contribution to be obtained, so other variables should 
be used. 

2 We assume that the non-resonant component is described by a smooth function. 
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Ks 




Fig. 3. The decay B ± -> K±t) c -> ^(KgR-nf. 




Fig. 4. The signal distribution of (A'si^vr 1 ' 1 ) invariant mass in the 
— )• K^^sKn) decay. The charmonium states r/ c , J/^>, x c i, and ij c (2S) (in 
order of mass) can be seen. The solid histogram is the combinatorial background 
determined from the AE sideband region. 

The Dalitz plots of 3-body r\ c and rj c (2S) decays are shown in Fig. 5. In the r] c 
case the distribution is not uniform and has a peaking structure around a Kit 
mass of 1.4 GeV/c 2 , which can be a combination of the iT*(1410), 1^(1430), 
and ^2(1430) states. However, the statistics are rather low, so it is impossible 
to determine with acceptable accuracy either the mass and width of these 
states, or the product branching fractions of r\ c decay modes involving these 
states. The small number of events does not allow the Dalitz analysis to be 
efficiently performed and makes it difficult to use the Dalitz plot variables to 
distinguish the r] c signal and non-resonant amplitudes. The same conclusion 
holds for the r] c (2S) decay. 

Another variable that can be used for the amplitude separation is cos 9. Since 
the T] c and r} c (2S) are pseudoscalars (J p = CT), we expect a uniform distri- 
bution in cos 9. In Figs. 6 and 7 the cos 9 distributions are shown for the rj c 
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Fig. 5. The Dalitz distributions in the r] c (left) and rj c (2S) (right) signal regions. 

and f] c (2S) signal and sideband regions, while the combinatorial background is 
subtracted. One can see that the sideband distribution has contributions from 
higher angular waves. A good fit can be obtained with the sum of S-, P-, and 
D-waves. The signal region also contains non-resonant background but mostly 
consists of signal events, so the S-wave contribution here is dominant. Separa- 
tion of the P- and D-waves from the S-wave in the non-resonant background 
allows the uncertainty from interference to be reduced. 




Fig. 6. Approximation of the cos 6 distribution in the rj c signal region by an S-wave 
(left) and in the sideband region by a sum of S-, P-, and D-waves (right). The 
combinatorial background is subtracted. 



Thus, we analyze a 2-D M(KKstt)-cos8 histogram assuming that the non- 
resonant signal amplitude is constant within the (2.5 - 3.46) and (3.14 - 
4.06) GeV/c 2 mass ranges. The number of events in a single bin is Nbi n = 
N/\e signal — k ■ N& Esideband , where the coefficient k is used to normalize the 
number of events in the AE sideband region. We minimize the likelihood 
function, assuming that the events in the signal and sideband regions are 
described by the Poisson statistics. In the r] c analysis the bin size along the 
cos# axis is 0.2 (9 bins), while along the M(KsK7r) axis it is 10 MeV/c 2 in the 
signal region and 150/130 MeV/c 2 in the left /right sideband region (44 bins). 
For the T) C (2S) the bin size along the cos 9 axis is also 0.2 (9 bins), while along 
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Fig. 7. Approximation of the cos 8 distribution in the rj c (2S) signal region by an 
S-wave (left) and in the sideband region by a sum of S-, P-, and D-waves (right). 
The combinatorial background is subtracted. 

the M(KsKtt) axis it is 16 MeV/c 2 in the signal region and 130 MeV/c 2 in 
the sideband region (29 bins). 

We exclude the J ft/) region (3.07 - 3.13 GeV/c 2 ) from the fit because the inter- 
ference of the J/ip and non- resonant signal is negligible due to the small width 
of the former and inclusion of this region does not constrain the i] c interference 
contribution. Moreover, the J/ip angular distribution has contributions from 
several amplitudes that are not well determined. The same arguments apply 
to the exclusion of the Xci mass region (3.48 - 3.54 GeV/c 2 ). We perform 
separate fits to the J/ip and xa with Gaussian functions using 1-D KgKn 
invariant mass distributions from signal MC and data. After comparing the 
obtained widths we determine the degradation of the resolution in data. Tak- 
ing this into account, we recalculate the detector resolution obtained from 
signal MC in the r] c and r) c (2S) regions. We obtain <j(r] c ) = (6.2 ± 1.1) MeV 
and a(r] c (2S)) = (9.8 ± 1.7) MeV. 

The fitting function can be represented as the square of the absolute value of 
the sum of the signal and non-resonant amplitudes integrated over all variables 
except M(K s Ktt) and cos 9: 



F{s,x) 



x 2 S 2 



N 



s' - M 2 + iMY 



1 i / i /2\ 

1 + E\X + E2X ) ■ 

M<ll <&) + a^sfe 2 , ql) J S{x') + 



PA P (q 2 1} 4)P{x') + 7 A D (& 4)D{x') 



ds'dx'dq l dq 2 d(j), (1) 



where x = cos 8, s = M 2 (KsKtt); q\ and q\ are Dalitz plot variables; E\ 
and £2 are constants that characterize the efficiency dependence on x and are 
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determined from MC; 5 and A are the bin widths in cos 6? and M(KsK7r) 
invariant mass, respectively; M and T are mass and width of the r] c (r] c (2S)) 
meson; N is the r] c (r] c (2S)) signal yield; a, /3, 7 are the relative fractions of the 
S-, P-, and D-waves, respectively; S = P = \f\x, D = fy^f (x 2 — |) are the 
functions characterizing the angular dependence of the S-, P-, and D-waves, 
respectively; A n is the signal S-wave amplitude, A StPjD are the background S-, 
P-, and D-wave amplitudes, respectively. The absolute values of the amplitudes 
squared are normalized to unity: 

J J J " \A v , s , P , D {ql,ql)\ 2 dqldqld4> = l. (2) 

To account for the momentum resolution, Eq. (1) is convolved with a Gaussian 
detector resolution function that is determined from the MC and calibrated 
from the J/tp (Xci) width in data. 

This function is determined by 15 parameters: N, M, T, a, /3, 7 described 
above; 6 parameters (JR v s, 9^s, $t v p, £%p, ^d, £%d) characterizing the inter- 
ference between the signal amplitude A n and background amplitudes As,p t D] 
3 parameters (Hsp, ^-sd, H-pd) describing the contributions from the interfer- 
ence between the background amplitudes As,p,d- In particular, 

$r,i + i%i = J J J A v (ql,ql)A*(qlql)dqldqld(t), (3) 

Ua = JJJ^ (A(ql ql)A*{ql q 2 2 )) dq\dq\d^ (4) 
where i, j — S,P,D and i ^ j. 

Since the function F(s, x) is a sum of the r\ c (rj c (2S)) Breit-Wigner and S-, P-, 
and D-waves, it can be represented as a rational function of s and x: 

In its most general form, such a function has 15 independent coefficients (Cy) 
in the numerator and two (M and V) in the denominator, however, in our case 
some coefficients are not independent: 

1. C 03 = M 2 (M 2 + T 2 )C 23 , 3. Co4 = M 2 (M 2 + T 2 )C 2A , 

2. Ci 3 = -2M 2 C 23 , 4. C14 = -2M 2 C 24 . 

Thus, we have (15 + 2 — 4 = 13) independent terms only, which is not enough 
to determine all 15 parameters of the function F(s,x). Two (15 — 13 = 2) of 
the parameters must be either obtained from other measurements or allowed 
to vary over the full allowed range. Since we have no additional information on 
these parameters, we scan over them. The result does not depend on the choice 
of the two scanned parameters. Since the interference is more significant in the 
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S-wave, we choose a and 5s v s- To perform this scan, we randomly sample a 
and Q v s in a reasonable rangdjD with the remaining 13 parameters free. After 
fitting the distributions we obtain a set of parameters and a x 2 • No additional 
local minima are found. 

The dependences of the signal yields on x 2 are shown in Fig. 8. In the rj c case, 
one can see that this distribution has a "plateau", which consists of fits with 
different N signa i (and other fit parameters) and the same x 2 ■ This feature arises 
because our system of equations for the fit parameters is underdetermined. 
The variation of parameters within this plateau will be referred to as the 
model uncertainty of our analysis and the average of their statistical errors 
as the statistical uncertainty. In the r] c (2S) case, the minimum x 2 plateau 
is not reached, because the parameters ^s^p^d, ^s^p^d, and Usp,sd,pd 
tend to their bounds[3- These bounds make our system of equations fully 
determined, so we can float the parameters a and 5s v s- In that case, the 
model and statistical errors cannot be separated. Thus we obtain N S i gna i = 
920±50(stat) ±170(model) for the ry c decay and N signa i = 128tll(stat+ model) 
for the T] C {2S) decay. 

WOO r 
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As A* /v niin 

Fig. 8. The dependences of the number of signal events on x 2 f° r the rj c (left) and 
rj c (2S) (right) decay analyses. One should note that the x 2 axis of the left plot 
covers a much smaller range than that of the right plot. 

Projections of the fits using the function F(s, x) are shown in Figs. 9 and 10. 

The fit procedure described above was also applied to MC signal samples. The 
obtained number of signal events was used to determine the detection efficiency 
((9.32 ± 0.10)% and (10.18 ± 0.10)% for rj c and rj c (2S) decays, respectively) 
and hence to calculate the product branching fractions. 



3 The scan range of parameter a is chosen so that it includes the minimum x 2 
region. It is [0:9] for the rj c and [1:11] for the r] c (2S) analyses. The parameter Q v s 
is varied over the entire physical region [-1:1]. 

4 By definition, $t v s, v P, v D, ^^s^p^d, and Usp,sd,pd vary in the interval [-1:1]. 
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2.5 2.75 3 3.25 3.5 -1 -0.5 0.5 1 -1 -0.5 0.5 1 



M(K s Kn), GeV/c 2 cos9 cos8 

Fig. 9. Projections of the fit in KsKir invariant mass in the r] c mass region (left) 
and cosf? in the rj c invariant mass signal (center) and sideband (right) regions. The 
combinatorial background is subtracted. The gap near 3.1 GeV/c 2 is due to the J/tp 
veto. The bin size along the cos 9 axis is 0.2. Along the M(KsKtt) axis the bin size 
is 10 MeV/c 2 in the signal region and 150/130 MeV/c 2 in the left/right sideband 
region 




M(K s Kn), GeV/c 2 cos9 cos9 

Fig. 10. Projections of the fit in KsKtt invariant mass in the rj c (2S) mass region 
(left) and cos 9 in the rj c (2S) invariant mass signal (center) and sideband (right) 
regions. The combinatorial background is subtracted. The gap near 3.5 GeV/c 2 is 
due to the Xci veto. The bin size along the cos 9 axis is 0.2. Along the M(KsKtt) 
axis the bin size is 16 MeV/c 2 in the signal region and 130 MeV/c 2 in the sideband 
region. 

The decay B ± — > K ± (KsKti) has two possible final states: K ± KsK ± n T and 
K ± KsK T 7i ± . We assume that the r\ c decay signal amplitudes are the same, 
but the non-resonant contributions (a, (3, 7) could be different in each de- 
cay channel. Because of limited statistics we do not treat these final states 
separately. A single distribution (1) effectively describes the incoherent sum 
of two distributions for the two decay channels. This does not affect the pa- 
rameters of r] c (i] c (2S)) states, but the parameters describing the non-resonant 
part obtained in the fit take effective values that depend on both non-resonant 
amplitudes. 

The method described above was checked using toy MC, which showed that the 
described procedure gives parameter values consistent with the generated ones. 
Moreover, a generic MC test was performed that included a full simulation of 
all b — > c decays at the T(4S) without interference. This test verified that the 
signal determination is not biased and gave an interference value consistent 
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with zero. 



4 Systematic uncertainties 

We evaluated possible sources of systematic uncertainties in the product branch- 
ing fractions. The number of BB pairs is calculated from the difference of the 
number of hadronic events on resonance and the scaled number of those off- 
resonance. The systematic error is dominated by the uncertainty in the scale 
factor and is equal to ~1.3%. We assume that the combinatorial background 
can be parameterized with a first-order polynomial. To obtain the background 
shape uncertainty, we describe the background by a second-order polynomial 
and compare the results. The uncertainty on the Kg decay branching fraction 
is taken from [6]. The contribution of the Ks reconstruction uncertainty was 
estimated in the Belle experiment to be 4.4% [21]. In our fitting procedure 
we take into account the efficiency dependence on cos 9 and assume that it 
does not depend on the KgK-n invariant mass. By adding a linear dependence 
on M(KsK7r) we estimate the corresponding systematic error. Moreover, we 
take into account the dependence of the efficiency on the r\ c and r) c (2S) de- 
cay models, such as KK*, KKq(1A30), and 7flf|(1430). The corresponding 
contribution to the systematic uncertainty is estimated by varying the effi- 
ciency obtained using these models and taking the difference in the results. 
An analysis of the charged track reconstruction uncertainty as a function of 
particle momenta has been performed in Belle data and gave an estimate of 
1% per charged track. To determine the errors due to K and ir meson identifi- 
cation, data from analysis of the process D* + — > D°tt + followed by the decay 
D° — > K~n + were used. The uncertainty in K ± identification is 0.8% per K 
meson and the corresponding value for 7r =1= identification is 0.5% per n meson. 
We also take into account the deviation of MC from the data by applying a 
correction to the efficiency: is 0.9996 for each kaon and 0.9756 for each 

J SMC 

pion. We vary the bin size along the x axis from 0.15 to 0.225, the AE window 
from 20 MeV to 40 MeV, and the detector resolution within the limits of its 
statistical uncertainty. 

Sources of systematic uncertainties for masses and widths include the back- 
ground parameterization, bin size, and detector resolution as described above, 
as well as a scale uncertainty and the effect of identical kaons in the final state. 
The scale uncertainty is determined from a comparison of masses of the J/ip 
(Xci) resonances, which were obtained by fitting the K s Kn invariant mass 
distribution, with the world average values [6]. In case two kaons of the same 
charge are present in the final state, the charmonium amplitude is a sum of 
two amplitudes corresponding to KsK^n and KsK^ir combinations. This 
can lead to the deformation of the K$K7r invariant mass distribution in the 
region of phase space where the KgK^n and K S K{ 2 )^ invariant mass values 
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overlap. We use toy MC to estimate this effect and take it into account as an 
additional systematic error. 



All the contributions to the systematic uncertainties are listed in Table 3 for 
the product branching fractions and in Table 4 for the masses and widths. 

Table 3 

Systematic uncertainties of the product branching fractions (in %). 



Source 


B ± -+ K- 






fir 


n c (2S) 


_ 

Number of BB pairs 


1.6 


1.3 




n 1 
U.l 


n 1 
U.l 


Model efficiency dependence 


+8.6 
-6.7 


+2.0 
-1.5 


JDdCKgl (J L1I1LL dJJJJI UJvlIIlfiLlUll 




4-9 ^ 


Bin size 


-3.3 


+ 13.3 
-3.9 


AE cut 


-2.2 


+2.3 


Detector resolution 


+1.1 


+4.7 
-8.6 


Mi nv efficiency dependence 


+2.2 


+0.8 


Track reconstruction 


3 


3 


identification 


1.6 


1.6 


-k^ identification 


1.5 


1.5 


Ks reconstruction 


4.4 


4.4 


Total, % 


+10.7 
-9.8 


+ 15.8 
-11.9 



Table 4 

Systematic uncertainties of masses and widths of the r] c and i] c (2S) mesons (in 
MeV/c 2 ). 



Source 


V 


c 


Vc(2S) 




Mass 


Width 


Mass 


Width 


Background approximation 






+0.2 


-0.1 


Bin size 


+0.2 


-1.0 


-1.1 


+2.4 


Detector resolution 


-0.1 


+ 1.0 
-1.2 


+0.5 
-0.1 


+ 1.0 
-0.9 


Scale uncertainty 


-2.0 




-1.7 




Effect of identical kaons 


+0.5 


+0.3 


+0.5 


+0.3 


Total, MeV/c 2 


+0.5 
-2.0 


+ 1.0 
-1.6 


+0.7 
-2.0 


+2.6 
-0.9 
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5 Results and discussion 



Table 5 shows a comparison of the results obtained assuming no interference 
(1-D fits to the AE and to the K s Kn invariant mass distributions) and those 
obtained using the analysis described above. One can see that taking interfer- 
ence into account leads to the introduction of a model error for the product 
branching fractions B(B ± — > K ± r] ( )B{ri c — > KsK ± tt t ) (for the i] c mass and 
width this error turns out to be negligibly small). In the r] c (2S) decay analysis 
the model error is not listed separately but the results differ noticeably from 
those that assume no interference. 



Table 5 

Comparison of the results obtained under the assumption of no interference between 
the signal and the non-resonant contribution and those obtained with interference. 





No interference 


Taking interference into account 


B ± K ± Vc: Vc (K s Ktt)° 


BxB, 1(T 6 
Mass, MeV/c 2 
Width, MeV/c 2 


24.0 ± \.2{stat)t 2 2 l{syst) 
2984.8 ±1.0{stat)±%l(syst) 
35A±3.6(stat)+ 3 2 ° 1 (syst) 


26.7 ± \A{stat)t 2 2 l{syst) ± A.9(model) 
2985.4 ± 1.5{stat)t°2 5 {syst) 
35.1±3.1(stat)ti° 6 (syst) 


B ± K±ri c (2S), Vc(2S) -> (KsKirf 


BxB, 1(T 6 
Mass, MeV/c 2 
Width, MeV/c 2 


3A±0.8(stat) ±0.2(syst) 
3646.5 ±3.7(stat)tH(syst) 
41.1 ± 12.0(stat)tio. 9 (syst) 


3.4+l 2 5 (stat + model)t i(syst) 
3636.ltH(stat + model^lisyst) 
6.6±H{stat + model)tH(syst) 



Table 2 shows that there is a large spread in the r) c (2S) width values. A 
possible explanation of this spread is that the previous studies did not take 
interference into account. For each of the studied processes the interference 
could have a different effect on the results and shift the r] c (2S) mass value 
significantly. Thus it is important to take interference into account. 

In addition to affecting the value of the branching fraction, interference changes 
the Breit-Wigner shape. This effect can allow the improvement of the statisti- 
cal accuracy with which the Breit-Wigner width is determined. In particular, 
the 7] C (2S) width, obtained in the present work, has a rather good accuracy, 
despite limited statistics and a detector resolution broader than the intrinsic 
width. The interference deforms the Breit-Wigner, lengthening its tail and 
thus improves the fit to the width (see Fig. 10). 
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6 Conclusion 



We report a study of the decay B ± — > K ± (cc), where the (cc) state decays 
to (KsKn) and includes the r\ c and r] c (2S) charmonia states. Both decay 
channels contain B ± — > K ± (KsKit) decays without intermediate charmonia 
that interfere with the signal. For the first time, the analysis takes interfer- 
ence into account with no assumptions on the phase or absolute value of the 
interference. 

As a result, we obtain an estimate of the model error for B(B ± — > K ± r] c )B(ri c — > 
K S KV): B(B ± ->■ K^Birjc -> K S KV) = (26.7 ± lA(stat)t 2 2 i(syst) ± 
A${model)) x lO" 6 . For B(B ± ->■ K^^S^B^^S) ->■ KsK^ir*), the model 
error from interference is not listed separately: B(B ± —> K ± rj c (2S))B(rj c (2S) — >■ 
K S KV) = (3Atii(stat + model)t° i(syst)) x 10~ 6 . 

We also obtain the masses and widths of i] c and r) c (2S). For the r\ c meson 
parameters the model error is negligibly small: 
M(t] c ) = 2985.4 ± l.5(stat)t°2 5 (syst) MeV/c 2 , 
r(r} c ) = 35.1 ± 3.1(stat)t\i(syst) MeV/c 2 . 

For the r) c (2S) meson the model and statistical uncertainties cannot be sepa- 
rated: 

M(r] c {2S)) = 3636. l±l%(stat + model)?%Z{sy st) MeV/c 2 , 
T(r] c (2S)) = 6.6±H(stat + model)±li(syst) MeV/c 2 . 

For the r) c (2S) the interference has a dramatic effect on the width (see Table 5). 

These results are consistent with those obtained in the most accurate existing 
measurements. Our errors are comparable to those in other experiments de- 
spite the fact that they include additional uncertainty related to interference 
effects. 
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Appendix 



After raising the absolute value to the second power F(s, x) (Eq. 1) convolved 
with the resolution Gaussian function can be written as: 



F(s, x) = £ / (1 + £l x' + e 2 x' 2 ) S 2 (x') (l m (s) + a 2 A + 2VNaI vS (s)) + 



2 



P 2 (x')/3 2 A + D 2 (x') 1 2 A + 2S(x')P(x')f3 (VNI vP (s) + aU SP A) + 



2S(x')D(x') 1 (yNI VD (s) + aU SD A) + 2P(x')D(x')/3 1 U PD A dj® 



where we use the following notations: 



2 



exp( 



'-w =///£? ? 



— M 2 + iMT 



A^qlq 2 ) ® 



V2na 



ds' ' dq\dq\d$ 
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1 



s' -M 2 + iMT 



6XP( ^ ds'dqldqldcP 



2na 



-f \s' -M 2 +iMT 



-f \s' — M 2 + iMT 



2na 



+ 



27TCT 



(8) 



I vP (s) 



3? 



s' — m 2 + /.wr 



^(g 2 ,g 2 2 )^(g 2 ,g 2 2 ) 



6XP( ^ ds'dqfdq^ 



c> c> 



V s' - M 2 + iMT ) s/2ia 

r ♦ ( \ 

Js-% \s' — M 2 + iMT 



+ 



(9) 



3? 



s' — m 2 + /.ur 



2tT(7 



ds'dq 1 dq 2 d(f) 



s+4 



27TfJ 



- M 2 + iMTJ 



! ) 

- M 2 + iMTJ 



ds' 



+ 



exp(-^r) 



27T(7 



CV2 



^5 = / / / H (A^i, <Z 2 2 )A* (g 2 , g 2 )) dg 2 dg 2 d0 = £„ 5 0T- 

^ S = / / / ^ Kfe 2 ^2 2 )^fe 2 ^2 2 )) dg 2 dg 2 d0, 
U v p = J J j U (A v (ql, qt)A* P (ql, qfj) dq\dq\d<\) = ^ P COS0„ P , 
%p= f J / 9 52)Ap(9i» 9 2 )) dg 2 dg 2 d0 = ^ P sin^ P , 

3?r,D = fff® (4,(9?, iD^Divh Ql)) dq\dq\d<\) = £ qD cos 0„ D , 



(10) 

(11) 
(12) 

(13) 

(14) 

(15) 
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%d = J J / 9 (A,(??> q 2 2 )A* D (ql ql)) dqldqldcj) = ^ D sin 9 vD , (16) 
UsP = I I j^{ A s(qlql)A* P {qlql))dqldqld<P, (17) 

UsD = I 1 1 ^{ A s(qlql)A* D (qlqi))dqldqld<p, (18) 

UpD = 1 1 1** ( Mq ^ & A M> d <& d <& d( t>- ( 19 ) 

Thus, the function F is determined by 15 parameters listed in Table 6. Pa- 
rameters a and $$ v s are fixed in the r\ c case, so they do not have a statistical 
error. In case the model error is not quoted, it is negligibly small compared to 
the statistical one. 

Table 6 



Fit results. 



Parameter 


T) ~\~ , r, ' -f- , / r, ' 7V" \ 

B ->■ if Jfc, Vc -> (K s Ktt) u 


B ->■ K*r) c (2S), r) c (2S) -)• (K s Ktt) u 


TV 


920 ± 50(stat) ± 170(modeZ) 


128±ll(stat + model) 


a 


4.8 ± 2.3(model) 


8At\^{stat + model) 


P 


5.3±1.6(stat)±2.1(model) 


3.7tH{stat + model) 


7 


7.9 ± 1.4(stat) 


6.2±\l(stat + model) 


£r?S 


0.37 ± 0.24(stat) ± 0.18(mocfeZ) 


0. 84tolt (stat + model) 




-0.06 ± 0.80(modeZ) 


-0.39tlil{stat + model) 


£r?P 


0.48 ± 0.22(stat) ± 0.19(mocfeZ) 


0. bit oil(stat + model) 


^P, rad 


0.08 ± 0.65(stat) 


-0.24tHt(stat + model) 


ir,D 


0.22 ± O.W(stat) ± 0.02(mode/) 


0. 8lt oll(stat + model) 


9 r)D , rad 


0.64±0.80(stat) 


0.58tiil(stat + model) 




-0.74 ± 0.46(stat) ± 0.23(mocfe/) 


— I^q (stat + model) 


n OT 


-0.11 ± 0.34(stat) ± 0.42(mocfe/) 


— I^q (stat + model) 


UPD 


0.24 ± 0.77 (stat) ± 0.50(mode/) 


-0.95+o^ (stat + model) 


M, MeV/c 2 


2985.4 ± 1.5(stat) 


3636.1^4^ (stat + model) 


T, MeV/c 2 


35.1 ± 3.1 {stat) 


6.6ttt(stat + model) 



Function F(s,x) can be represented as a rational function of s and x: 



„2 2 4 



(20) 



where 
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(1) Om = f 7 2 , 

(2) C 23 = IVTEfaiipD, 

(3) C 22 = §/3 2 - f 7 2 + iVbarflLsD, 

(4) C 21 = v^a/OT^ - 

(5) C 20 = \a 2 + | 7 2 - ^« 7 n 5 D, 

(6) C14 = -2M 2 C 24 , 

(7) Ci 3 = -2M 2 C 23 , 

(8) Ci 2 = IVEVN-ytonD ~ 2M 2 C 22 , 

(9) C n = y/Sy/NptonP - 2M 2 C 21 , 

(10) Cio = v^VaS^s - ^fVN^ D - 2M 2 C 20 , 

(11) C 04 = M 2 (M 2 + r 2 )c 24 , 

(12) C 3 = M 2 (M 2 + r 2 )C 23 , 

(13) C 02 = |>/5v^V : 7(Mr^ r , - M% D ) + M 2 (M 2 + r 2 )C 22 , 

(14) Cm = ^^(Mrv - M% P ) + M 2 (M 2 + r 2 )c 21 , 

(15) C 00 = |JV + y/Na(MT% s - M 2 ^ s ) - ^VN^{MF% d - M 2 % lD ) + 

M 2 (M 2 + r 2 )c 20 . 
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